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ABSTRACT 



Context. Subdwarf B pulsators exhibit pressure (p) and/or gravity (g) modes. Their frequency spectra range from very simple, with 
few frequencies, to very rich, with more than fifty peaks in some cases. Balloon09 is a hybrid pulsating subdwarf B, showing a large 
number of p- and g-modes including a triplet and a quintuplet-like structures which are interpreted as p-mode frequency splittings 
due to stellar rotation. Photometric observations undertaken in two subsequent years revealed a change in these rotational splittings 
of 0.24- 0.58 pHz/yr. 

Aims. We analyse the possibility of angular momentum interchange between stellar rotation and internal gravity waves as a mechanism 
for the observed rotational splitting variations. 

Methods. The expected change in the rotational splitting of eigenmodes resulting from this mechanism are computed in the non- 
adiabatic linear approximation for a stellar structure model that is representative of the target star. The fact that g-modes are also 
observed in Balloon09 makes it a particularly suitable candidate for our study, because the change in the rotational splittings are 
proportional to the amplitude squared of the g-modes which, in this case, can be estimated from the observations. 
Results. We find that this mechanism is able to predict changes in the splittings that can be of the same order of magnitude as the 
observed variations. 

Key words, subdwarfs - stars: oscillations - stars: rotation 



1. Introduction 

Subdwarf B (sdB) stars populate the blue extension of t he hor- 
izonta l branch (EHB) in a Hertzsprung-Russell diagram dHeberi 
I1986I) . Their high temperatures (r eff > 20 000 K) and high grav- 
ities (5 < log g < 6) correspond to an evolutionary state in be- 
tween the red giant and white dwarf phases, although the de- 
tailed history is not well understood yet. An anomalously high 
mass loss rate at some point is needed to explain the hot sub- 
dwarf's thin H-envelope, wh ich prevents them from ascending 
the as ymptotic giant branch ( iDorman et al. Cruz et alj 

Il996l) . Mass transfer due to binary interaction has been in- 
voked as an explanation, as it is observationally supported by 
the high binary fraction of these objects in the field. The pre- 
dom inantly single sdBs that populate the EHB in globular clus- 
ters dMoni Bidin et al.ll201 ll) can be explained by white dwarf 
mergers (Han et al. 2008). 

Some sdB stars exhibit stellar oscillations. Indeed, long- 
(~lh) and short-period (~ 1 min ) pulsating sdBs ca n be dis- 
tinguished dKilkennv et all 119971: iGreen et al] l2003[). with a 
few hybrid objects show i ng both period r e gimes JSchuh et alj 
l2006t lOreiroetalJ 12005 iLutz et alj l2009t iBaran et al.l I20TTT 
Long-period oscillations are ex plained in terms of grav- 
ity modes dFontaine et al.l |2()()3). wh ile short-period varia - 
tions are attributed to pressure modes dCharpinet et alj|1997l) . 
Asteroseismic techniques can thus be applied to sdBs to re- 



trieve information on their internal structure and rotation, and 
will even tually help to constrain their evolutionary formation 
channels dHu et alj|2008l) . 

Balloon 090100001 (Bal09 hereafter) is a hybrid sdB pul- 
sator, with a rich frequency spectrum both at low and high 
frequencies. Together with the facts that it is the brightest 
sdB pulsator and has one of the highest amplitude of oscil- 
lation, make Bal09 a very inte resting object. Sinc e the dis- 
covery of its pulsating nature dOreiro et al. I I2001 . a multi- 
site photometric campaign was organized to unravel the com- 
plex frequency structur e of the target during the 2005 sum- 
mer (IBaran et al.l 120091). while a low resolution time-resolved 
spectroscopy ( Telting & 0stensenl 120061) and a hi gh resolution 
spectroscopy dataset en abled a mode identification dTelting et alj 
l2008tlBaran et al.l2008l) . Mode id entification was also att empted 
through multicolour pho tometry dCharpinet et alj|2008l) . while 
Ivan Grootel et a fl d2008h provide a seismic solution for Bal09. 

An interesting feature concerning Bal09's amplitude spec- 
trum is that a triplet-like structure is resolved close to the dom- 
inant, single mode. Moreover, five components are also iden- 
tified ~25 //Hz apart from the triplet. From these modes a fre- 
quency splitting of ~ 1.5 /iHz is derived. Assuming a canonical 
total mass M = O.5M and using the spectroscopic log g - 5.39 
value , a rotational velocity between 1.5 - 2km/s dTelting et al.l 
l2008h is obtained, in agreem ent with the typic al low surface- 
rotation speed of single sdBs dGeier et al ] l2010h . 
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Bal09 has been photometrically monitored after the main 
multi-site campaign in 2005 but no result has been published 
yet (Baran, private communication). A comparison of the 2004 
and 2005 seasons, allowed the detection of changes in the ampli- 
tudes of many modes. Amplitude variatio ns are not uncom mon 
among pulsating sdBs, as it is discussed in lKilkenrivl(l2010l) . On 
the other hand, small changes in the frequency splitting are found 
from one campaign to the other as well. Although close to the 
frequency resolution, the triplet components have moved in the 
right direction as if rotation is varying. The quintuplet also varies 
in the same way, despite the fact that not all five components 
were resolved in 2004. Since it is unusual for a star to change 
its rotation in such a short time scale, a physical explanation for 
this phenomenon is yet to be found. 

The aim of this work is to examine one possible mechanism 
that could explain the changes in the frequency splitting. We 
analyse here the possibility of angular momentum interchange 
between stellar rotation and internal waves. This scenario re- 
quires the existence of internal gravity waves with a dissipation 
mechanism or any other non-conservative process that eventu- 
ally could lead to amplitude variations. This is the case in Bal09, 
where we can estimate the g-mode energies from the observa- 
tions and hence quantify the angular momentum transport by 
waves. 

The possibility of angular momentum interch ange b e tween 
wave s and rotation was originally proposed by lAndol (119811 
119831) . where it is described how non-axisymmetric, non-radial 
oscillations can redistribute angular momentum in stars. The 
complete formalism is included in these works, and is also ap- 
plied to two Ceph stars. He concluded that non-radial os- 
cillations can affect significantly the rotation profile of these 
stars, but in a time scale (10 4 yr) much larger than the one ob- 
serve d in our case. A related process has been postulated for the 
Sun (IZahn et al.lli997t iKumar et al.lll999t iTalon etalj|2002b to 
explain the fiat Sun's rotation profile in the radiative interior. In 
that case however, the angular momentum is carried out by very 
low frequency gravity waves generated at the base of the con- 
vection zone, which are then completely damped in the radiative 
interior. The time scale f or such a moment um transfer is about 
10 7 yr for solar-like stars dTalon et al.ll2002l) . 



2. Observational data 

In Fig. Q] we compare the amplitude spectrum of the target ob- 
tained from the 2004 and 2005 photometric campaigns. The 
figure shows the residual spectra after the three highest ampli- 
tude p-modes (at ~2.8 mHz) have been prewhitened. In general, 
lower amplitudes are obtained for p-modes in 2005, whereas 
similar or larger are obtained for g-modes. This suggests that 
p- and g-modes are interchanging energy, a fact that could be 
related to the changes in the rotational splittings, but in any case 
the analysis of this process is outside the scope of the present 
work. The highest amplitude g-modes for the two photometric 
runs are compared in Table [TJ A complete compa rison of modes 
in the two sessions is found in Ba ran et al. 1 (120091) . 

In Fig.|2]a closer view to the the triplet and quintuplet struc- 
tures is shown as observed in 2004 and 2005. Interestingly, 
these changes do not seem random, but as if the mean rota- 
tional velocity of the star was increasing in the time scale of 
the observations. The retrograde \m\ = 1 triplet component 
changed by 0.24 pHz/yr from 2004 to 2005, while the retrograde 
\m\ = 2 quintuplet component changed its frequency as much as 
0.58 pHz/yr in the same period. 
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Fig. 1. Comparison of the amplitude spectrum of Bal09 in 2004 
and 2005 after r emoving the 3 highest peaks. Figure from 
iBaran et alj d2009l) . P, g, and c stand for p-modes, g-modes, and 
combination frequencies respectively. 

Table 1. Frequency analysis for the highest amplitu de g-modes 
based on photometric data. From lBaran et al. 1 J20091) . 



2004 


2005 


Freq. 


Ampl. 


Freq. 


Ampl. 


(pHz) 


(mma) 


G"Hz) 


(mma) 


229.57 


0.57 


229.55 


0.85 


239.97 


1.96 


239.97 


1.08 


246.31 


0.61 


246.30 


1.52 


272.38 


2.75 


272.46 


1.72 


298.89 


0.67 






325.67 


1.25 


325.61 


4.44 


331.21 


0.81 


331.18 


0.75 


365.81 


2.66 


365.81 


2.65 


397.19 


0.49 


397.23 


0.84 


631.08 


0.44 


630.74 


0.42 


684.35 


0.60 


684.40 


0.47 


833.08 


0.53 


833.09 


0.43 



Telting & 0stensen (2006) publish frequencies and line-of- 
sight velocities based on low resolution mode spectroscopy ac- 
quired during the summer of 2004. Their results are included in 
Table [2] High-re solution spec t roscop y during the summer 2006 
was obtained bv lTelting et al.l (120081) . They only list the line-of- 
sight velocity for the fundamental mode: 14.5 km/s, which is sig- 
nificantly lower than 18.9km/s measured in 2004. In both cases, 
line-of-sight velocities were computed from cross-correlation 
profiles (ccp), which merge all the spectral lines into a sin- 
gle, and with higher signal-to-noise artificial line. Thus, spec- 
tral features formed at different depths are combined to pro- 
duce the ccp's. More over, B aimer lines a r e used in the low- 
resolution dataset by ITelting & 0stensenl d2006l) . while only 
m etallic ones are cros s-correlated from the high-resolution data 
of ITelting et al.l d2008l) . Hence, in this case, the change in ampli- 
tude can be attributed to the different techniques used. 

3. Structure models and non-adiabatic oscillations 

3.1. Model 1 

For the most part of the work we have considered a single stellar 
structure model that closely matches the fundamental parameters 
of Bal09 (Model 1 hereafter). It is a 0.43M Q total mass model 
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Fig. 2. Splittings of the € = 1, v = 2807.5 ,uHz (left panel) and 
t = 2 v = 2853 4 iMz (right panel) modes in 2004 and 2005. 
From lBaran et al.l (120091) . 

Table 2. Frequency analysis based on low-resolution spec- 
troscopy by Telting & 0stensen (2006). The observations were 
carried out in 2004. 



Frequency 


RV 


0"Hz) 


(km/s) 


272.44 


0.86 


325.67 


0.96 


365.77 


1.05 


2807.47 


18.89 


2823.24 


5.88 


2824.81 


3.44 


2826.28 


1.86 



with a hydrogen amount of 2 x 10~ 4 M G , r eff = 26 800 K, and 
log g = 5 .47. Note that the stellar mass equals the seismic mass 
found by Ivan Grootel et al] d2008i) . The T s ff and \ogg are within 
the observed spectroscopic errors considered by these authors, 
although we chose T e g on the low side to ensure unstable g- 
modes. 

The structure model is c onstructed with a version of the stel- 
lar evolution code STARS (Egg letonlll97lh . updated for aster- 
oseismology. We started the evolution at the ZAEHB, and the 
particular model we consider is near the end of core He burn- 
ing with an EHB age of 2.1 x 10 s yr. H/He diffusion is in- 
cluded consistently during the evolution, whereas the effects of 
diffusion of Fe and Ni are approximated with a Gaussian accu- 
mulation around log T = 5.3. Radiative opacities that account 
for the Fe/Ni-enhan cement are taken from the Opacity Project 
(iBadnell et al.ll2005l) . More details and input physics are given 
in Hu et al.l ( 20091) and references therein. 

Normal modes of oscillation and eigenfunctions for 6 < 6 
were computed in the line ar, non-adiaba tic approximation with 
the pulsation code MAD (iDupretl 1*200 ll) . The structure model 
was chosen so that the fundamental mode (v = 2.786 mHz) 
is close to Balloon's observed highest amplitude mode at v = 
2.807 mHz. This frequency does not show a ny rotational split- 
ting and has been identified as a radial m ode (iBaran et alj|2008t 
iTelting et al.ll2008t ICharpinet et al.ll2008l 

In Fig. [3] the Brunt- Vaisala' (N) and Lamb (Si) frequencies 
for this model are represented. They are displayed in terms of 
cyclic frequencies so that the propagation cavities are easily de- 
limited in comparison with observed modes. Peaks in N denote 
gradients in the composition profile, hence, we can distinguish 
the H envelope (from r/R ~ 0.45 to the surface), the He ra- 
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Fig. 3. The black solid line is the Brunt- Vaisala frequency, N, 
while red, green, blue, and yellow lines are the Lamb acoustic 
frequencies, S e, for t = 1, 2, 4, and 5, respectively. The dashed 
horizontal lines correspond to the modes indicated in the text 
(the I = 1 and { = 2 p-modes look like a single thick line). 
Model 1 was used. 

diative layer (0.2 < r/R < 0.45) and the inner convective CO 
core formed as ashes of the He-core burning process. The hori- 
zontal lines in the figure correspond to the frequencies of some 
selected modes: two /7-modes at v = 2807.5 yt/Hz and v = 2853.4 
//Hz, which are the 6=1 and ( = 2 theoretical counterparts of 
the observed splittings, and two low frequency g-modes: t = 4, 
y = 436.0 //Hz and I = 5, v = 29 1 .4 //Hz. 

In the left panel of Fig.|4]we show all the computed eigenfre- 
quencies for this structure model, indicating if the mode is lin- 
early unstable (in green) or not (in black). Low order p-modes 
(2800 < v < 4000 //Hz) are excited for all the angular degrees {() 
considered, which can account for the high-frequency peaks de- 
tected in Bal09. Another island of unstable modes is due to high 
order t > 3 gravity modes. Their frequencies, within 400 - 900 
//Hz are slightly shifted compared to the observed low-frequency 
peaks in Bal09. The relevance of this fact in our work will be 
discussed later on. 

The radial and horizontal components of the displacement 
eigenfunction vector for some particular modes (those indicated 
in Fig. [3} are included in Fig. [5] The contribution of every point 
in the star to the dimensionless kinetic energy, &, given by 



An^pr 2 [t(€+im\ 2 + \f r \ 2 ]dr 



M[((e+m h (R)\ 2 + \UR)\ 2 } ' (1) 

is also included in Fig [5] The first two modes on the top are 
the theoretical eigenfrequencies representative of the observed 
triplet and quintuplet respectively. The I = 4 and the t = 5 are in 
the observed g-mode region, albeit only the I = 4 is theoretically 
unstable. 

3.2. Model 2 

One of the problems with the current sdB models is that the 
observed instability strip of the g-modes can not be correctly 
modeled. Compared to the observations, typically, the T^'s of 
models with unstable g-modes are too low, the periods of un- 
stable g-modes are too low and only mo des with high spher- 
ical degree (t > 3) are excited (s ee e.g iFontaine et al.|l2003t 
Uefferv & Saiol2006tlHu"et al.l2009l) . This problem is also appar- 
ent in our Model 1 as discussed in the previous section. However, 
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Fig. 4. Theoretical eigenfrequencies as function of angular de- 
gree. Modes with angular degrees < ( < 6 and cyclic frequen- 
cies v > 150 //Hz were considered. Excited modes are shown in 
green. Left panel is for Model 1 and right panel for Model 2. 
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Fig. 5. Radial (red, dotted line) and horizontal (green, dot-dashed 
line) components of the eigenfunctions normalized such that 
% r {R) = 1 (p-modes) or %h(R) = 1 (g-modes), where R is the 
photospheric radius, for some oscillation modes: From top to 
bottom: i = 1, v = 2807.5 //Hz, £ = 2, v = 2853.4//Hz, € = 4, 
v = 436.0 //Hz and I = 5, v = 29 1.4 //Hz. Also shown (black 
solid line) is the integrand in the dimensionless kinetic energy & 
defined by Eq.Q] Model 1 was used. 



Hu et al. (in prep.) have computed new models with atomic dif- 
fusion, including radiative levitation of H, He, C, N, O, Ne, Mg, 
Fe and Ni, that solve the aforementioned issues with the unsta- 
ble g-modes. Hence, we also compare with such a preliminary 
model (Model 2 hereafter). Model 2 is obtained from the same 
ZAEHB model as used for Model 1 but now the evolution and 
pulsations are computed with fully self-consistent atomic diffu- 
sion. 



The right panel of Fig |4]shows the stable and unstable modes 
for Model 2. In this case there are a large number of unstable 
modes in the observed frequency range, between 200 //Hz and 
800 //Hz and including degrees from t — 1 to 6. Note that for 
higher degrees the frequency range of unstable modes is shifted 
upwards. 

4. Transfer of angular momentum and change in the 
rotational splitting 

Non-radial, non-adiabatic oscillations, especially g-modes with 
large tangential velocities, can transfer angular momentum to the 
mean rotational flow. This variation in the rotational velocity will 
consequently affect the frequency splitting of non-radial g and 
/?-modes. 

4. 1 . Transfer of angular momentum due to gravity modes 

The angular momentum transfer due to nonaxisymmetric waves 
in sta rs can be described in terms of a Reynolds stress dAndol 
I1981I) . namely: 



-(pm 2 £l) 
ot 



-t where t = V C ■ (pmVpV^,) 



(2) 



Here m, V p , V$ and V; are the cylindrical polar radius, poloidal 
component of the perturbed velocity, azimuthal component of 
the perturbed velocity, and poloidal component of the differential 
operator V, respectively. The rotation angular velocity is given 
by Q.. Here and in the following we will use a similar notation to 
lUnno et al.l(ll989h . 

Taking the linear and Cowling approximations in the wave 
equations, assuming a dependence of the form exp [i{crt + m<p)], 
and if the star rotates slowly, such that the real part of the eigen- 
frequencies obey <r R » f2, the change in the angular momentum 
can be expressed as: 



T — r wave + TNA 

where the terms T wave and tna are given by: 



= cr R crimp 



tna = -mP\m 



dPYIdp 
PJ\P 



/in 1 2 



(3) 



(4) 



(5) 



Here cr = <x R + i<T\ is the complex eigenfrequency of the mode 
and Yf the spherical harmonic. Tna represents the transfer of 
angular momentum due to non-adiabatic effects. In this approx- 
imation tna is equal to the work function, w(r), except for a 
constant factor. On the other h and, the term T wave represents a 
wave transience. As explained in lAndoldl98ll) . in Eq.|?]for T wave , 
the imaginary part of the eigenfrequency, err, comes from a time 
derivative of the wave amplitude and hence can be derived ap- 
propriately. In particular, it makes sense to take <ti from the ob- 
served changes in the amplitudes of the g-modes, provided they 
are intrinsic and not caused by beating phenomena. 

Fig. [6] shows r wave and tna as given by Eq. |4] and [5] respec- 
tively for an t — 4 unstable mode. According to Eq. [2] these 
functions give the contribution of each point in the star to the 
change in the angular momentum. The wave transit term, T wave , 
is similarly weighted as the integrand in the dimensionless en- 
ergy S (compare to third panel in Fig. |5j, while tna, even for 
inner gravity modes, is only important in the uppermost layers 
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Fig. 6. T wave and tna normalized to the maximum value for a 
I = 4 g-mode with v = 436.0 //Hz. For Tna only the external 
layers are shown since it is negligible at deeper ones. Model 1 
has been used. 

where non- adiabatic effects are large. The relative importance 
between both terms depends on the energy of the mode and the 
value of <Ti considered in r wave and will be discussed in the next 
section. 

Since sdB stars are slow rotating stars, in principle Eqs. [4] 
and [5] would be enough for estimating the changes in rotation. 
However, prograde and retrograde modes contribute in an oppo- 
site way to the angular momentum transfer. If modes with the 
same degree but azimuthal orders +m had the same amplitude, 
the net effect would cancel out. For this reason we will also con- 
sider the full expression for tna as given by Eq. 36.16 in Unno 
etal. (1989): 



T NA f = pmlm 
— mpRe 



(cr + mLl)\'T — 

C n 



(v p -p- l V t p)v T 5S*lc p 



(cr* + mQ) 



(6) 



where v p , v$ and 6S are linear perturbations in, respectively, the 
poloidal and azimuthal component of the velocity, and the en- 
tropy. This equation does not neglect £1 against cr R nor does it 
assume the Cowling approximation. However, we shall still con- 
sider the eigenfunctions of a non-rotating equilibrium model. 

4.2. Effect in the rotational splitting 

The rotational splitting, i.e. the pe rturbation in the f requencies 
caused by rotation, is given by (e.g. lAerts et al.ll2010L Eq. 3.349) 



6&nim -ml I K n ( m {r, #)Q(r, 6)drd8 
Jo Jo 

where the kernel K„( m (r, 6) is given by; 
sm0{\t r (r)\ 2 P>;Xco S 6) 2 



(7) 



K n tm — 



-Mr)Y 



dP'l'Y m 2 

1 + _^p;»(cos<9) 2 

sin" 
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Fig. 7. Contribution of the different layers to the change in the 
splitting of the I = 1, v = 2807.5 //Hz mode due to the t = 4, 
m = 4, v = 436 /iHz g-mode according to Eq. [TO] Top panel is 
for T wave and bottom panel for tna- Values have been normalized 
to the maximum. Note that only the upper layers of the model 
are shown since the contribution of deeper layers is negligible. 
Model 1 has been used. 



-P?{cosff) 2 [£{r)t h {r) + f,(r)^(r)] 



-2PT(cosey 



dP" cosO. 



d6 sin6» 



pr /I„lm 



and 



Inlm — 



2 (7+M)! 

2zTT(z- \m\y. 



\trf + l(l+ 1)161' 



pr dr . 



(8) 



(9) 



Therefore, from Eqs.|2]and|2l we obtain the variation of the fre- 
quency splitting caused by transfer of angular momentum be- 
tween eigenmodes and rotation: 



<9w 



'nCm 



dl 



Jr*R r-n 
I K n ( m 
o Jo 



*n'l'm' 

pr 2 sin 2 6 



-drdO 



(10) 



where «, I, and m are for the mode for which the rotational split- 
ting is considered and «', C , m' are for the mode for which the 
transfer of angular momentum is computed. 

Fig. [7] shows the normalized outcome of Eq. \W\ when ei- 
ther the wave transit term T wav e (top panel) or the non-adiabatic 
term tna (bottom panel) are considered. Here we have chosen 
the splitting of the theoretical t — 1 p-mode closest in frequency 
to the observed triplet and computed the change due to a single, 
unstable, t = 4 g-mode. For T wave , the change in the splitting is 
mainly caused in the external layers, whereas the angular mo- 
mentum transfer takes place in deeper regions, see Fig. [6] This 
is a straightforward consequence of the nature of the /j-mode 
considered for measuring the rotational splitting that have larger 
amplitudes in the external layers. For tna, the external layers 
are dominant for both the rotational splitting variation and the 
angular momentum transfer, although the weights are different. 
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4.3. Visibility 

Since all computations are based on the linear theory, the eigen- 
functions are calculated except for a constant factor that can be 
fixed by, for example, the velocity amplitude at a given point 
in the atmosphere. In principle velocity observations are eas- 
ier, less model dependent, and more related to theoretical values 
than photometric observations, and therefore we will use them. 
However, as commented in Sect. 2, the line-of-sight-velocity ob- 
servations do not correspond to a specific point in the atmo- 
sphere. Since the amplitude of the modes change by more than 
one order of magnitude throughout the atmosphere, an accurate 
determination would require a detailed theoretical simulation of 
the observational measurements. In this work we will perform a 
simple calibration between observed and theoretical amplitudes 
by assuming that the observations correspond to a unique optical 
depth, tr, but considering different values. 

Specifically, we use a relation between the amplitude of the 
radial component of the velocity, v p , and the observed line-of- 
sight velocity, v b s . The proportionality factor, that we will call 
the visibility, A{£, m, w), is given by 

v bs = A(£, m, cJ) v p at the observed point. (11) 

Given that we observe v Q b s for certain modes, we can re- 
trieve their v p if we are able to compute A((,m,a>). The factor 
A({, m, cJ), includes the geometrical effects of the velocity pro- 
jection as well as limb darkening effects. It also depends on the 
inclination angle i and the ratio between the horizontal and verti- 
cal components of the eigenfunctions, K = at the observed 
point. To avoid the dependence on i, the rms value over m can be 
used. When referring to this average, we will just remove the m 
dependence, A(£, at). 

To compute ACL m, a >) we have mostly followed the formal- 
ism in lAerts et al.l(l 19921) . However, in that work the approxima- 
tion 

K ~ (GM)/(flV) (12) 

was used. This approximation is derived by assuming a 
Lagrangian pressure perturbation at the outer boundary of 6P = 
0, which is valid for adiabatic oscillations, provided that the 
boundary condition is not located too deep in the atmosphere. 
However, Eq. [T2]is not valid if non-adiabatic effects are impor- 
tant. Since we have computed the theoretical eigenfunctions in 
a non-adiabatic approximation for our specific models, we can 
compute the precise ratio K at the required point in the atmo- 
sphere. 

In Fig. [8] we show the K ratio evaluated at two optical depths. 
The approximate K from Eq. Q~2] is also included for compari- 
son. Note that the approximation is in better agreement with the 
actual theoretical values at higher depths, where non-adiabatic 
effects are less important. The visibility factors A(t, oS) depend 
on the optical depth through limb-darkening coefficients and the 
ratio K which is negligible for /9-modes. However for g-modes 
it can be as large as 100 or 1000 (as can be seen in Fig. [8]) and 
as a result g-modes can have rather different visibility factors 
depending on the optical depth considered. 

Since the horizontal components of the g-modes are very 
large, the visibility factors as defined by Eq.QTJare not by them- 
selves indicative of the kind of modes we expect to observe. It 
is more illustrative to express the visibility in terms of the en- 
ergy of the modes. The kinetic energy, av eraged over time, can 
be expressed as (see e.g. lAerts" et al.ll2010l) : 

£k,n = jBMVL = -&Mv 2 p {R) 1 1 + j (13) 
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Fig. 8. K ratio between horizontal and radial components of the 
displacement vector at optical depths log tr = -3 and log tr = 
-4. Black points include all the linearly stable modes while 
green points are for the unstable ones (Model 1). Red crosses 
are the ratios expected by the approximate Eq. [12] 

where & is given by Eq.[TJ V rms is the root mean squared velocity 
over the stellar surface, at radius R and M is the total mass of 
the star. Using Eqs. [T3l and [TT1 we obtain the following relation 
between the line-of-sight velocity and the kinetic energy of the 
modes: 

= A\l,m,u» 2 
E kin M8 1+K 2 £(£+1)' 

The right hand side of Eq.[l4]can be theoretically computed, 
and it is shown in Fig. [9] for two different optical depths. They 
have been normalized to the fundamental mode. Note that g- 
modes in the observed frequency range have much larger values 
than the fundamental one and hence need much less excitation 
energy to be observed with similar Doppler velocities. This is 
true not only for the lowest degrees but also for I = 4 and I — 5. 
Moreover this conclusion is valid for any optical depth in the 
range considered. 

4.4. Amplitude calibration 

According to Eqs. [3] |4j |5] and [10] the transfer of angular mo- 
mentum and the change in the splittings are proportional to the 
eigenfrequency squares and hence the amplitude square of the 
g-modes considered. Since the amplitude of the modes can in- 
crease by an order of magnitude throughout the atmosphere, a 
very important issue concerning Eq.[TT]is where to fix the point 
where obs erved veloci t ies are transformed to theoretical ones. 
Following iRauch et"aD (1201 Ot) . the formation depth of the core 
of the Balmer lines in the subdwarf stars seems to be at optical 
depths between tr ^0.1 and tr > 0.01. However the wings of 
the lines are formed higher in the atmosphere and since the am- 
plitude of the modes increase with height, the oscillatory signal 
could come from optical depths substantially lower than that of 
the core. Here we shall try not to overestimate the change in the 
angular momentum due to the oscillations. This requires to as- 
sume as small as possible g-mode energies, and hence to locate 
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Fig. 9. Ratio |v b s |/ V^kin normalized to the fundamental mode. 
See text and Eq.[14]for an explanation. Modes with the same de- 
gree have been joined for clarity. Black, red, green, blue, yellow 
and magenta correspond to angular degrees £ — 0, 1, 2, 3, 4, and 
5, respectively. Upper panel is for log r R = -3 and bottom panel 
is for log t r = -4. 

the signal as high as possible in the atmosphere. We will take 
tr = 0.001 as a reference value. 

Since only few g-modes are observed in velocity and they 
cannot be identified (/ and m are unknown), we need some hy- 
pothesis about the excitation energy to be able to proceed. To 
calibrate the theoretical amplitudes we attribute a line-of-sight 
velocity of v Q b s = 1 km/s to the g-mode with the largest v? b /Skin 
in the frequency range [250, 400] ji/Hz. The value of v Q b s and the 
frequency range considered are based on the observations. Our 
choice gives the minimum energy required to reproduce the ob- 
servational amplitudes in the g-mode region. At the two optical 
depths considered in Fig. [8] the theoretical mode that matches 
this requirement is a £ — 3, v = 276 //Hz eigenmode for Model 
1. but there are other I = 2, ... ,5 modes that if taken as repre- 
sentative of the observations would result in a similar E^\ n , so 
the particular mode is not important for an order of magnitude 
estimation. 

On the other hand, the optical depth of the model at which 
the line-of-sight velocity is calibrated has a large impact on the 
results. For the above g-mode, the amplitude of the radial veloc- 
ity eigenfunction at an optical depth of log tr = -3 is 5.3 times 
smaller than that at log tr = -4. Since the change in the splitting 
is proportional to the velocity square, a factor of 28 difference is 
found between both optical depths. 

5. Results 

In the following computations we assume the same kinetic en- 
ergy for all modes, although there is no particular theoretical rea- 
son for this. However, from the velocity observations, we known 
that at least a few modes in the frequency range [250, 400] fiHz 
need to have such a large E^ n and the observations in photometry 
reveal that in fact a bunch of g-modes are excited to similar am- 
plitudes and hence, perhaps, energies. On the other hand, if £kin 
varies in orders of magnitude with frequency or angular degree, 
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Fig. 10. Change in the rotational splitting of the £ = 1, v = 
2807.5 yuHz p-mode due to the £ — 4, m — 4 g-modes, consid- 
ering the tna term. Here we assume the same energy for all the 
modes and fix it so that v b s = 1 km/s for £ = 3, v = 276 ^/Hz at 
log tr = -3. Black points are for stable modes and green points 
for the unstable ones. Model 1 was used. 



the actual behaviour of the changes in the rotational splittings, 
proportional to E\^ a , can be very different from those shown in 
the following figures. Hence, we can estimate the order of mag- 
nitude of the change from individual modes but we are not able 
to compute a net effect from all the g-modes. 



5.1. The non-adiabatic term, t na 

Let us start by considering the changes in the splitting due to 
the non-adiabatic term tna of individual g-modes with m > 0. 
In Fig. [10] we show the change in the splitting of the £ = 1, 
v = 2807 //Hz mode due to g-modes with £ = 4, m = 4. Eq. [5] 
for Tjva and Model 1 have been used. In the figure, the energy 
£kin was fixed by taking the value v b s = 1 km/s at log tr = 
-3 for the reference mode given in Sect. 14.41 The modes that 
induce the largest changes in the splitting are in the observed 
frequency range, although most of them are not expected to be 
excited in Model 1. It follows from Fig.[l0]that the change in the 
splitting of the theoretical counterpart of the observed triplet is as 
large as 0.02 /iHz/yr if log tr = -3 is assumed, considering only 
the contribution of the £ = 4, m = 4 g-modes. Slightly larger 
values are obtained if £ = 5, m = 5 g-modes are considered. 
These changes for single modes are only an order of magnitude 
smaller than the observed values. On the other hand, as stated 
before, calibrating the observed amplitude at an optical depth 
of log tr = -4 results in changes in the rotational splittings of 
about 28 times smaller. Thus the resulting splitting for the same 
single modes would be about 0.007 ^Hz/yr. 

However the net effect could be much smaller since modes 
with the same degree £ and radial order n but different azimuthal 
orders m can be expected to be excited to very similar amplitudes 
and, according to Eq. |5j the change in the rotation due to ±m 
pairs would mostly cancel out. Note also that high frequency 
g-modes give a contribution opposite to that of low frequency 
g-modes, as seen in Fig.fTOl 

To examine the possible reduction of the net effect, we now 
assume the same amplitude for all the 1i + 1 azimuthal modes 
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Fig. 11. Change in the rotational splitting of the I = 1, v = 
2807.5 //Hz p-mode due to ( = 4 multiplets, assuming the 
same amplitude for all the modes with different m but the same 
I, n numbers, and considering the TNA f term. As before we as- 
sume the same energy for all the modes and fix it so that the 
I = 3,v = 276 //Hz mode would have an observed velocity 
v b s = 1 km/s if it were to be observed at an optical depth 
log tr = -3. Model 1 has been used. 



with the same angular degree I and radial order n and use the full 
expression of TNAf (Eq. |6j to compute the change in splitting. 
Fig. fTTI shows the changes in the splitting of the same t = 1, 
v = 2807.5 //Hz p-mode due to the € = 4, g-modes. The energy 
of the modes were again fixed by assuming an optical depth of 
log tr = -3. In this case the change induced by isolated (2€ + 1) 
group of modes can be as large as 0.004 //Hz/year, which is 60 
times smaller than the observed one. Again, considering the ( = 
5 modes, slightly larger values are obtained while for log tr = 
-4 the effect will be reduced by a factor 28. 

We could sum up the contribution of all the modes thus prob- 
ably making the change in rotational splitting even larger than 
we give here. However, as explained at the beginning of Sect. 5, 
we do not find it realistic to give a conclusive summed up num- 
ber. Instead, our results should be considered as a low order of 
magnitude estimation of the change in rotational splitting due to 
angular momentum transfer. 

5.2. The transient term, T wave 

Here we consider the changes due to the term T wave . This term 
is caused by a change in the mode amplitude, represented in Eq. 
[4] by the term cr t . Since the eigenfrequencies are computed in 
the linear theory, considering the theoretical values cr\ has little 
sense. On the other hand, as shown in Table Q] the mode ampli- 
tude is changing yearly between campaigns. These changes can 
be caused by different factors, such as beating phenomena, but as 
an upper limit we can assume they are completely due to intrinsic 
changes in the mode amplitudes. A value of cr\ — 1 /year is con- 
sistent with the results in Table [T] In any case, since the change 
in the splittings are proportional to cr\ the results are easly trans- 
formed to other amplitude variations. Figure[T2]shows the results 
in the splitting of the 6 = 1, v = 2807.5 //Hz p-mode due to the 
t — 4, m — 4 g-modes. Compared to Fig. [10] we see that this 
changes are about 100 times smaller than those corresponding 



Fig. 12. Change in the rotational splitting of the I — 1, v = 
2807.5 //Hz p-mode due to the I — 4, m = 4 g-modes consid- 
ering the T wave term. Here we assume the same energy for all the 
modes and fix it so that v bs = 1 km/s for I = 3, v = 276 //Hz at 
log tr = -3. Black points are for stable modes and green points 
for the unstable ones. Model 1 has been used. 



to Tna- Thus, unless the changes in the mode amplitudes have 
time scales much shorter than a year, the contribution of T wave 
can be neglected as compared to the tna effect. 

5.3. Changes in the splitting for different modes 

So far we have limited the computations to the changes in the ob- 
served triplet (assuming a ( = 1, low order p-mode). However, 
the change in the splitting strongly depends on the frequency of 
the p or g-mode used as the test. On the other hand, the expected 
change for the 1 = 2 modes are very similar to that of the I = 1 
when compared at a given frequency, except for a factor m for 
p-modes and a factor m{\ - + 1)]} for g-modes, both well 
known results from the asymptotic theory. To illustrate the fre- 
quency dependence, Fig. Qj] shows the changes in all the I = 1 
modes due to a g-mode with I = 4, v ~ 300 //Hz. From Fig. 
[T3l follows that the largest changes are expected for p-modes 
with high radial order. The observed triplet is around 2800 //Hz 
but higher p-modes are in fact observed. So this result is poten- 
tially very interesting since we predict changes in the splitting 
around or above 1 //Hz/year. Also, the change in the rotational 
splitting of low frequency g-modes would be significantly larger 
than those corresponding to the observed modes, but measuring 
rotational splitting in this region could be a difficult task due to 
the large population of peaks and beating phenomena in the sig- 
nal. 

Since tna is proportional to the work function, the change in 
the splittings are determined by the same non adiabatic proper- 
ties of the g-mode eigenfunctions as the growth rates, although 
the integrals over the whole star have different weights. Since 
the modes with the largest contribution to the changes in the ro- 
tational splittings are stable in Model 1, this can be an important 
source of errors. On the other hand, as shown in Fig. [4] g-modes 
with / < 5 and low frequencies, are unstable in Model 2 and 
hence, in principle, using this new model would avoid that prob- 
lem. Fig.[T3l allows to analyze this problem. The g-mode chosen 
for the transfer of angular momentum is expected to be stable 
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Fig. 13. Change in the rotational splittings of the 6 = l,m = 1 
modes due to the and ( = 4, m = 4, v ^ 300 jt/Hz mode. Absolute 
values have been considered. The energy of the mode was fixed 
in such a way that v b s = 1 km/s for ( — 3, v — 276 //Hz at 
log tr = -3 for Model 1. Black point are for Model 1 and red 
points for Model 2. The g-mode considered has v = 305.14ji/Hz 
in Model 1 and v = 301.02 yt/Hz in Model 2. 



if Model 1 is used and unstable if Model 2 is used. Despite this 
fact, Fig. [13] shows that the results for both models are rather 
similar, or a little larger for Model 2, if compared at the same 
frequency of the test mode. This allows us to be confident in the 
results obtained with Model 1 . The reason for using Model 1 in 
the previous figures was that the frequencies of the / = 1 and 
/ = 2 lowest /7-modes are very close to the observed triplet and 
quintuplet while, as previously mentioned, we do not yet have a 
model with the new physics that accurately match the observed 
p-modes. 

6. Conclusions 

Transfer of angular momentum between g-modes and rotation 
can explain the changes in the rotational splittings observed in 
the sdB star Balloon 090100001 . We have done the computations 
by estimating the mode energies from the observed g-mode ve- 
locities and using structure models that fit the stellar parameters, 
including the frequencies of the highest amplitude p-modes. The 
theoretical changes found for the rotational splittings due to sin- 
gle g-modes are about one order of magnitude smaller than the 
observed ones if an optical depth of log tr = -3 is considered 
for the calibration. If the cancellation effect between prograde 
and retrograde modes are considered by assuming the same am- 
plitude for the rotational multiplets, the net effect is about 60 
times smaller than the observed ones. Since the optical depth 
considered is rather a lower limit and because the observed value 
should be compared with the net effect from all g-modes, we can 
conclude that this mechanism can produce a rotational splitting 
of the same order of magnitude as the observed ones. Moreover, 
we find that the transfer of angular momentum is possible due 
to the non-adiabatic nature of high radial order, low degree g- 
modes, even if the amplitude of the modes do not vary. In fact, 
changes in the amplitudes within a yearly time scale, as it seems 
to be the case in this and other sdB stars, induce a much smaller 
change in the rotational splittings than the non-adiabatic process. 



Although Ba09 has brought us an opportunity for checking 
this kind of stellar phenomena, some uncertainties prevent us 
for obtaining more quantitative results. One uncertainty comes 
from performing an accurate calibration between observed ve- 
locities and theoretical energies, because the former cannot be 
easily translated to values at an specific optical depth or to a 
some weighted integral over the atmosphere. Photometric obser- 
vations can help here, but a calibration of such observations are 
generally more complex and model dependent than the velocity 
ones since changes in the flux are due to non-adiabatic effects. 
This problem could be limited in the future by performing de- 
tailed theoretical simulations of the data analysis. 

A second problem that hinders more quantitative calcula- 
tions concerns the mode identification in the low frequency 
range. But even if this were possible the main contribution to 
the transfer of angular momentum can come from high degree 
modes that are not observed due to geometrical cancellations. 
However it is interesting to note that potentially observed modes 
and modes contributing the most to the changes in the rotational 
splitting can be rather close in frequency and degree. As can be 
seen in Fig. [4] Model 2 predicts unstable modes with a frequency 
range that resemble the observed one. There are unstable modes 
for all the angular degrees considered, but for higher degrees the 
unstable modes shift towards higher frequencies. Thus the main 
contribution, that according to Figs. [10] and [TT] are for frequen- 
cies below 400 or 500 //Hz could come from modes of degrees 
in the observed range or at most a little higher, say t < 5 or 6. 

According to Eq. [2] since r depends on latitude, so does 
the change on the rotational velocity. Hence, g-modes induce 
a differential rotation in the stars. It is interesting to note that 
the quintuplet did show a latitudinally differential rotation, and 
hence the same angular momentum transfer can be invoked as 
the cause for the differential rotation. 

Searching for other sdB stars with similar characteristics 
would allow to improve our knowledge. This in volves in par- 
ticula r space mission like COROT and KEPLER dBorucki et al.l 
1201 Ol) that can acquire high precision and long-term photome- 
try of pulsating hot subdwarfs, which is of utmost importance to 
detect the possible effects in frequency splittings due to angular 
momentum transfer. However, t here is only one p-m ode pulsator 
identified in the Fo V of Kepler dKawaler et alj20fob and a hand- 
ful of g-mode pulsating sdBs with candidate s plitting patterns 
that s t ill nee d longer time series to be confirmed (I0stensen et al.l 
l20loil20TTl) . 
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